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To  explore  the  dynamic  characteristics  of  the  SOFC  systems  and  to  develop  relevant  control  strategies,  a 
previously  developed  steady  state  SOFC  model  is  converted  to  a  dynamic  model.  The  model  includes  mass, 
momentum,  thermal  and  electrochemical  analysis,  as  well  as  the  kinetic  model  of  hydrocarbon  reactions. 
Applying  two  control  strategies  i.e.,  cell  constant  fuel  flow  rate  and  constant  fuel  utilization  during  the 
transient  time,  the  model  is  implemented  to  analyse  the  dynamic  behaviour  of  a  planar  direct  internal 
reforming  (DIR)  SOFC  cell  under  several  step-load  changes.  Transient  response,  resulting  from  an  inlet 
temperature  variation,  is  also  investigated.  The  results  show  that  the  relaxation  time  is  strongly  related 
to  the  thermal  behaviour  of  the  cell  and  the  applied  control  strategy.  However,  it  is  almost  independent 
of  the  load  variation  magnitude. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

With  the  advantages  of  its  high  efficiency,  high  reaction  kinet¬ 
ics,  high  power  density,  fuel  flexibility  and,  as  well  as  high  level 
of  usable  heat;  the  high  temperature  solid  oxide  fuel  cell  (SOFC) 
is  identified  as  one  of  the  most  promising  technologies  for  future 
power  generation  in  both  large  central  power  station  plants  and 
distributed  power  generation  stations.  However,  the  SOFC  tech¬ 
nology  faces  many  challenges  when  it  comes  to  commercialization, 
since  cost  reductions,  increasing  the  reliability  of  components  and 
extended  durability  are  required.  To  overcome  these  challenges, 
one  needs  to  understand  the  complex  operating  mechanism  of 
the  SOFC  stack,  such  as  heat  and  mass  transfer,  together  with 
electrochemical  reactions.  Under  steady-state  operation,  the  high 
temperature  and  aggressive  environment  inside  the  SOFC  stack,  as 
well  as  the  failure  of  fuel  cell  system  components,  can  cause  harmful 
conditions  on  the  SOFC  stack.  In  addition  to  steady  state  operation, 
as  the  demand  for  electrical  power  is  frequently  fluctuating,  the 
technology  also  has  to  remain  reliable  and  stable  with  transient 
situations  such  as  start-up,  shut-down  and  load  change.  Therefore, 
it  is  increasingly  important  to  understand  the  fundamental  mech¬ 
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anisms  and  SOFC  behaviours  under  transient  operating  conditions. 
Furthermore,  knowledge  about  the  SOFC  dynamic  behaviour  can  be 
employed  to  design  a  control  scheme  for  the  SOFC  operation  and  to 
compensate  its  dynamic  interactions  with  the  rest  of  the  system. 

As  the  experimental  investigations  on  SOFCs  are  time- 
consuming  and  expensive,  mathematical  modelling  techniques 
have  been  proved  to  be  a  cost  effective  method  to  investigate 
both  steady  state  and  transient  behaviour  of  SOFCs.  There  are 
many  examples  of  modelling  investigations  in  literature.  Most 
models,  however,  are  in  steady  state  condition  [1-5]  and  tran¬ 
sient  SOFC  models,  typically,  are  not  as  elaborate  as  steady-state 
models.  Achenbach  presented  one  of  the  first  transient  studies  on 
SOFCs  [6].  He  considered  a  three-dimensional  and  time-dependent 
mathematical  model  of  a  planar  SOFC,  together  with  semi  empir¬ 
ical  voltage  polarization  equations.  Ota  et  al.  [7]  investigated  the 
dynamic  behaviour  of  a  micro-tube  SOFC.  A  dynamic  model  for 
tubular  SOFC  developed  by  Sedghisigarchi  and  Feliachi  [8]  based 
on  the  lumped  capacitance  model.  Aguiar  et  al.  [9]  presented  the 
open  and  closed-loop  transient  response  of  a  co-flow  planar  anode- 
supported  intermediate-temperature  (IT)  direct  internal  reforming 
(DIR)  solid  oxide  fuel  cell  to  load  step  changes.  A  reduced  ID 
dynamic  model  of  a  planar  DIR-SOFC,  with  assumption  of  uniform 
current  density  distribution  and  one  temperature  layer,  for  system 
research  is  introduced  by  Kang  et  al.  [10].  A  thorough  review  of  the 
state-of-the-art  in  dynamic  modelling  is  presented  in  [11,12]. 

The  overall  aim  of  the  present  study  is  to  evaluate  the  perfor¬ 
mance  of  a  planar  SOFC  system  for  residential  applications  in  both 
steady  state  and  transient  conditions.  Results  of  the  cell  and  system 
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Nomenclature 

Latin  letters 
A  area(m2) 

AER  excess  air  flow  ratio 

Q  molar  concentration  of  component  i  (mol  m-3 ) 

CP  specific  heat  (kj  kg-1  K-1 ) 

%ernst  Nernst  cell  potential  (V) 

E  rate  of  energy  flow  (kj  s-1 ) 

^Proud,  i  rate  °f  energy  transfer  accompanying  mass  transfer 
of  product  i  (kj  s-1 ) 

Ereact,  i  rate  of  energy  transfer  accompanying  mass  transfer 
of  reactant  i  (kjs-1) 

FU  fuel  utilization  (%) 

H  hight  (m) 

i  cell  current  (A) 

i'o  exchange  current  (A) 

J  current  density  (A  m-2 ) 

LHV  lower  heating  value  (J  mol-1 ) 

h  molar  flow  rate  (mol  s_1) 

hf  consumed  fuel  consumed  in  each  channel  (mol  s-1 ) 

PDC  stack  DC  power  output  (W) 

q  heat  flux  (kW  m-2 ) 

fj  rate  of  reaction  j  (mol  m-2  s-1 ) 

T  temperature  (I<) 

t  time  (s) 

Vop  cell  voltage  (V) 

V  volume  of  a  computational  element  (m3 ) 

x  axial  coordinate  (m) 

Greek  symbols 

rjact  activation  polarization  (V) 

?7cen  cell-stack  efficiency 

diffusion  polarization  (V) 

Pohm  ohmic  polarization  (V) 

Vij  stoichiometric  coefficient  of  component  i  in  reaction 

j 

p  density  (kg  m-3) 

r  relaxation  time  (s) 

Superscripts 
cond  conduction 

conv  convection 

Rad  radiation 

Subscripts 
A  air 

An  anode 

ave  average 

Ca  cathode 

El  electrolyte 

F  fuel 

Final  final  condition  after  variable  change  (new  steady 
state  condition) 

Initial  initial  condition  before  variable  change  (old  steady 
state  condition) 

IcA  air  side  interconnector 

IcF  fuel  side  interconnector 

In  inlet 

Nernst  Nernst  potential 

ohm  ohmic 

out  outlet 

Ox  oxidation  reaction 


PEN  positive-electrode/electrolyte/negative  electrode 

red  reduction  reaction 

SR  steam  reforming  reaction 

WGS  water  gas  shift  reaction 


models  in  steady  state  condition  have  been  previously  published 
in  [13,14].  In  this  article,  as  an  incremental  step  toward  a  more 
comprehensive  understanding  of  the  dynamic  behaviour  of  planar 
DIR-SOFC  systems,  a  previously  developed  steady  state  SOFC  cell 
model  is  converted  to  a  dynamic  model.  By  using  this  model,  the 
SOFC  dynamic  behaviour  is  investigated  based  on  two  strategies 

1. e.,  constant  flow  rates  to  the  system  and  constant  fuel  utilization 
during  the  cell  transient  operation.  Transient  response,  resulting 
from  an  inlet  temperature  variation,  is  also  investigated. 

2.  SOFC  cell  and  stack  model  description 

Only  a  brief  description  of  the  cell  models  applied  is  given  here 
(for  details  see  [13,14]). 

The  considered  SOFC  configuration  is  a  planar  type  in 
both  anode-supported  (or  intermediate  temperature  (IT))  and 
electrolyte-supported  (or  high  temperature  (HT))  designs  (Fig.  1). 
A  general  detailed  quasi  2D  model  based  on  the  mass,  momen¬ 
tum,  thermal  and  complete  electrochemical  analysis,  as  well  as  the 
kinetic  model  of  hydrocarbon  reactions,  is  developed  for  the  con¬ 
sidered  geometry.  The  planar  SOFC  cell  is  divided  into  a  series  of 
piped-type  volumes  (unit  elements)  and  the  finite  volume  method 
is  used  to  discretize  the  governing  equations  on  each  unit  element. 
The  main  assumptions  of  the  cell  model  are  [13]: 

(1)  One-dimensional  cell  representation  along  the  streamwise 
direction. 

(2)  Uniform  distribution  of  feed  gases  to  each  channel. 

(3)  Each  of  the  gas  channels  in  the  unit  cell  act  as  continuously 
stirred  tank  reactors  (CSTR). 

(4)  Lumped  temperature  of  the  solid  cell  structure 
(anode-electrolyte-cathode). 

(5)  Adiabatic  boundaries  at  the  cell  inlet  and  outlet. 


Fig.  1.  Planar  SOFC  geometry. 
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Due  to  the  high  electrical  conduction,  electrodes  and  current 
collector  act  as  isopotential  surfaces. 

The  electrochemical  analysis  comprises  the  complete  evaluation 
of  the  cell  ohmic,  activation  and  concentration  polarizations.  These 
polarizations  are  subtracted  from  the  Nernst  potential  to  calculate 
the  unit  element  operational  voltage: 

V0p  —  ^Nernest  —  (^ohm  +  ^act,  Ca  +  *7act,  An  +  ^7 diff,  Ca  +  ^diff,  An)  0) 

For  calculation  of  the  ohmic  polarization,  equivalent  ohmic 
resistance  is  calculated  for  each  unit  element  by  dividing  it  into  sev¬ 
eral  parts.  Ohm’s  law  is  then  applied  to  calculate  the  ohmic  losses. 
As  the  reactant  gaseous  concentrations  in  the  triple  phase  boundary 
(TPB)  are  needed  to  calculate  the  concentration  polarization,  two 
levels  of  diffusion  phenomena,  between  the  bulk  gas  phase  to  the 
electrode  surface,  and  between  electrode  surface  and  triple  phase 
boundary  are  considered.  The  Fick’s  model  is  then  applied  to  cal¬ 
culate  the  reactant  gaseous  concentrations  in  the  electrode  surface 
and  TPB,  respectively.  The  Fuller  equation  is  used  to  calculate  the 
binary  diffusion  coefficients.  The  Knudsen  diffusion  coefficients  are 
calculated  based  on  the  average  molecular  speed  of  each  compo¬ 
nent.  The  effective  diffusion  coefficients  are  also  calculated  based 
on  the  binary  and  Knudsen  diffusion  coefficients  with  consideration 
to  the  mixing  rule. 

The  complete  form  of  the  Butler-Volmer  equation,  with  selec¬ 
tive  transfer  coefficient  (two  step  equation  [15]),  is  also  applied  to 
calculate  the  activation  polarization. 

The  continuity  equations  comprise  the  anode  and  cathode  molar 
flow  rate  variations  due  to  the,  internal  reforming  of  methane,  and 
water-gas-shift  reactions: 


=  qfcT^cond  -  qfc0A,Vconv  -  f/pEN,  icA^Rad  Air  side  interconnect 

(7) 

PlcF,  ave^IcFCPicF,  ave  ^ 

=  qfc°FndAcond  -  qfc°FnFAcond  -  dpEN,  icF^Rad  Fuel  side  interconnect 

(8) 

In  the  model,  the  fully  developed  laminar  flow  solutions  of  the 
Navier-Stokes  equations  or  steady  Hagen-Poiseuille  analysis,  is 
used  for  calculation  of  the  pressure  losses  inside  the  air  and  fuel 
channels  [16].  The  enthalpy  of  the  mixtures  in  the  flow  channels 
are  evaluated  by  the  molar  mixing  law.  The  kinetic  theory  of  gases 
is  used  to  calculate  the  dynamic  viscosity  and  thermal  conductivity 
of  the  mixtures. 

For  discretization  of  the  spatial  and  temporal  terms  in  the 
applied  equations,  the  implicit  upwind  and  second  order  central 
scheme  are  utilized,  respectively. 

3.  Model  parameters 

SOFC  fuel  utilization  (FU),  air  excess  ratio  (AER)  and  efficiency 
respectively  are  defined  as: 

yp  _  ^f,  consumed  /qn 

(4hCH4  +  nco  +  Hh2  )An>  In 


d(Q) 

dt 


~^  +  E  Vj  ''S<CH4’C°2< 

j= WGS,  SR,  Ox 


CO,  H20,  H2,  N2}  anode 


(2) 


+  5>rj  i e  {02,  N2)  cathode  (3) 

j=red 

As  described  in  [13],  to  improve  the  accuracy  of  the  thermal 
model,  five  layers  of  temperatures  are  considered  in  energy  bal¬ 
ance  equations.  Radiation  between  solid  parts  is  also  considered 
in  the  thermal  model.  Therefore,  for  each  unit  element,  the  energy 
balance  equations  are  written  for  the  fuel  and  air  side,  PEN  (positive 
electrode-electrolyte-negative  electrode),  and  air  and  fuel  streams 
as: 

PF,  ave^F^F,  ave  —  ^in,  F  —  ^out,  F  +  ^pen^F^001™  +  ^IcF^fcA^cond 

-  fteact,  h2  +  E proud,  h2o  Fuel  flow  (4) 


PA,  ave  VaCPa, 


dT\ 
ave  dt 


£  £  i  ^conv  a  i  ^cond  a 

Mn,  A  -  ^out,  A  +  Wpen,  A7^01™  +  ^IcF,  IcA^cond 

-%eact,o2  Airflow  (5) 


PPEN,  aveWENCPpEN,  ave 


dlpEN 

dt 


_  ncond a  ,  nconv  A 

—  4  pen  ^cond  WPEN, 


conv 


AER  = 


_ (rc°2  )ca,  in _ 

[2hCH4  +  0.5(hco  +  Hh2  )]An  In 


dCell  = 


Ppc 

(hp  LHVf)An  jn 


(10) 

(11) 


Relaxation  time  (r)  is  also  an  important  parameter  in  investigat¬ 
ing  the  dynamic  behaviour  of  fuel  cells.  It  is  defined  as  the  period 
necessary  to  recover  a  percentage  (usually  90%,  95%  or  even  99%) 
of  the  dynamic  voltage  drop  (or  increase)  resulting  from  the  dif¬ 
ference  A V=V{t^  oc)  -  V(t  =  +0)  [6],  where  the  t  =  +0  term  shows  a 
very  short  time  after  the  load  change. 


4.  Results 


This  section  is  divided  into  two  sub-sections.  Section  4.1  gives 
results  of  the  cell  transient  responses  to  several  load  changes,  for 
different  control  strategies  (first  strategy:  example  Cases  1  and  2, 
second  strategy:  Case  3).  Section  4.2  gives  results  of  the  cell  tran¬ 
sient  responses  to  an  inlet  flow  temperature  change  (Case  4). 

A  co-flow  planar  electrolyte-supported  DIR-SOFC  is  selected  for 
the  present  study.  The  input  parameters  of  the  model  are  shown  in 
Table  1.  As  can  be  seen,  methane  with  30%  pre-reforming  is  con¬ 
sidered  as  input  fuel  composition  to  the  model.  During  the  load 
and  temperature  changes,  the  AER  parameter  is  fixed  to  a  constant 
value.  It  does  not  necessarily  mean  that  the  air  flow  rate  is  com¬ 
pletely  constant.  Based  on  Eq.  (10),  in  a  constant  AER,  the  air  flow 
rate  is  dependent  on  the  fuel  flow  rate  to  the  cell. 

4.1.  Transient  response  to  step-load  changes 


^PEN^A^conv  +  9 pen,  IcA^Rad  +  QpEN,  IcF^Rad  +  ^React,  H2 
T^React,  02  —  ^Produ,  H20  —  PEN 


PlcA,  ave^IcA^PlcA, 


dT\CA 


ave 


dt 


In  this  study,  the  transient  behaviour  of  a  planar  DIR-SOFC  is 
investigated  under  two  control  strategies.  In  the  first  strategy,  it  is 
assumed  that  the  fuel  and  air  flow  rates  are  constant  after  the  aver¬ 
age  current  density  step-change  (or  step-load  change).  As  a  result, 
the  FU  parameter  changes  to  a  new  value  which  is  proportional 
to  the  new  average  current  density.  Therefore,  in  a  positive  load 


Table  1 

SOFC  model  input  parameters. 
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Geometry 
Anode  thickness/m 
Cathode  thickness/m 
Electrolyte  thickness/m 
Interconnector  thickness/m 
Cell  active  area 
(width  x  height)/mm2 
Channel  height,  fuel  side/m 
Channel  height,  air  side/m 
Channel  width,  fuel  side/m 
Channel  width,  air  side/m 
Number  of  channels,  fuel  side 
Number  of  channels,  air  side 
Material  electrical  properties 
Conductivity  of 
anode/ W  irr 1  K_1 
Conductivity  of 
cathode/Wnrr1 K-1 
Conductivity  of 
interconnector/W  irr 1  K-1 
Interconnector  specific 
resistivity/^2  m 

Cathode  specific  resistivity/^  m 

Electrolyte  specific 
resistivity/^2  m 

Anode  specific  resistivity/^2  m 
Operating  conditions 
Fuel  utilization/% 

Average  current  density /A  irr 2 
Air  excess  ratio 
Air  feed 
Fuel  feed 

Cells  outlet  pressure/Pa 
Fuel  inlet  temperature/0 C 
Air  inlet  temperature/°C 
PEN  emissivity 
Interconnector  emissivity 
Methane  reforming  data 
K  coefficient 
m  coefficient 
n  coefficient 
Methane  activation 
energy/J  mol-1 
Activation  polarization 
Activation  energy  of 
anode/kj  mol-1 
Activation  energy  of 
cathode/kj  mol-1 
Pre-exponential  coefficient  for 
anode/A  m-2 

Pre-exponential  coefficient  for 
cathode/A  m-2 
m  coefficient 
Diffusion  polarization 
Pore  diameter  of  anode/m 
Pore  diameter  of  cathode/m 
Porosity  of  cathode  /% 

Porosity  of  anode  /% 

Tortuosity  of  cathode 
Tortuosity  of  cathode 
WGS  reaction 
KWGS  (mol  s_1  m-2  Pa-1 ) 


50 x 10“6 
50 x 1CH5 
150 x 10-6 
2.5  xlO-3 
100 x 100 

1  x  10“3 
1  x  10“3 
3  x  10-3 
3  x  10-3 
18 
18 

2.0 

2.0 

2.0 

^oiexp(^) 

42.0x10®  (  —1200  ^ 

7>EN  eXP  V  TPEN  ) 

3-34xl04exp(=lSf5) 


85 

3000 

7 

21%  02,  79%  N2 

Steam/carbon  >  2, 30%  pre-reforming 

101,325 

900 

900 

0.8 

0.3 

4274 

1.0 

0.0 

82,000 


100 

117 

5.5  xlO”8 
7  x  10-8 
1.0 

1  x  10“6 

1  x  10“6 

50 

50 

3.0 

3.0 

0.01 


change  C/finai  -./initial >  0)*  the  cell  starts  with  a  given  low  fuel  utiliza¬ 
tion  and  reaches  a  higher  value  for  the  new  steady  state  operation. 
In  a  negative  load  change  C/finai  - ./initial  >0)>  the  cell  starts  with  a 
given  high  fuel  utilization  and  reaches  a  lower  value  for  the  new 
steady  state  operation. 

With  regard  to  the  first  strategy,  an  example  of  the  cell  tran¬ 
sient  response  is  presented  in  Fig.  2  (Case  1 ).  For  this  example,  the 
initial  average  current  density,  fuel  utilization,  cell  efficiency,  and 
operating  voltage  (at  t<0),  respectively  are,  5000  Am-2,  FU  =  85%, 
ijceii =  =39.7%  and  V0p  =  0.520  V.  At  t  >  0,  the  average  current  density 


Time ,  s 


28 

26 

24  * 
22  | 
o 

20 

18 

16 

14 


Fig.  2.  Transient  responses  of  cell  voltage  and  power  to  step-load  change  from 
5000 Am-2  to  3000 Am-2  (Case  1). 


Fig.  3.  Current  density  distribution  along  cell  length  as  a  function  of  time  (Case  1 ). 


is  changed  to  3000  Am-2.  As  can  be  seen  in  Fig.  2,  there  is  an  over¬ 
shooting  of  the  cell  voltage  above  the  new  steady  state  voltage.  This 
overshoot  has  been  observed  by  several  authors  [6,9,17-20]  and  it 
is  correlated  to  the  slow  thermal  response  of  the  SOFC  materials. 
The  overshoot  magnitude  reported  in  the  literature,  however,  is  not 
the  same  and  it  is  highly  dependent  on  the  applied  control  strategy, 
as  well  as  the  cell  variable  parameters. 

The  present  model  output  results  for  rjCen,  FU,  V(t=+0)  and 
V(t—>  oo )  are  28.3%,  51%,  0.687  and  0.618,  respectively.  The 


Fig.  4.  PEN  temperature  distribution  along  cell  length  as  a  function  of  time  (Case  1 ). 


8952 


P.  Kazempoor  et  al.  /  Journal  of  Power  Sources  196(2011)  8948-8954 


1350 


1300 


1250 


1200 

1150 

1100 

1050 


Fractional  cell  length 


200 


Vn 


400 


3 

2 


Fig.  5.  Fuel  flow  temperature  distribution  along  cell  length  as  a  function  of  time 
(Case  1). 


relaxation  times  for  recovering  90%,  95%  and  99%  of  the  dynamic 
voltage  increase  respectively  are  r90%  =  270s,  r95%  =  315s  and 
r99%  =  390s. 

Figs.  3,  4  and  5  show  the  current  density,  the  PEN  (positive 
electrode-electrolyte-negative  electrode)  and  fuel  temperature 
distributions  along  the  cell  length,  after  the  average  current  den¬ 
sity  step-change  (or  step-load  change),  until  the  new  steady-state 
operation  point.  Due  to  the  lower  average  current  density  and  the 
lower  FU  in  the  SOFC  new  operation  point,  a  smoother  distribu¬ 
tion  for  all  three  quantities  can  be  observed  in  the  new  steady  state 
condition.  Flowever,  in  spite  of  the  cell  voltage  variation  (Fig.  2),  no 
jump  in  the  PEN  and  fuel  temperature  distributions  is  observed. 

In  Fig.  6,  the  maximum  and  minimum  temperatures  of  the  cell, 
as  well  as  the  maximum  cell  temperature  gradient,  are  shown. 
As  can  be  seen,  the  maximum  allowable  local  temperature  gra¬ 
dient  and  the  maximum  allowable  total  temperature  difference 
( A^total  =  Tmax  —  ^min)  along  the  cell  are  lower  in  the  new  steady 
state  operation.  However,  high  variation  of  the  local  temperature 
gradient  and  ATtotal  may  cause  damaging  thermal  stress  within  the 
SOFC  components  Therefore,  although  the  above  described  strat¬ 
egy  is  an  easy  way  to  control  an  SOFC  system,  it  is  not  an  effective 
method  and  can  cause  harmful  conditions  during  the  SOFC  opera¬ 
tion  points. 

As  a  second  example  (Case  2),  Fig.  7  shows  the  cell  tran¬ 
sient  response  to  a  new  current  density  step-change  i.e., 
4000-3000 Am-2.  For  this  example,  the  model  output  results  for 
rj ceii,  FU,  V{t  =  +  0)  and  V{t^o c)  are  28.3%,  51%,  0.64  and  0.618, 
respectively.  The  relaxation  times  for  recovering  90%,  95%  and 
99%  of  the  dynamic  voltage  increase  respectively  are  r9o%  =  300s, 
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Fig.  7.  Transient  responses  of  cell  voltage  and  power  to  step-load  change  from 
4000 Am-2  to  3000 Am-2  (Case  2). 

r95%  =  350  s  and  r99%  =  430  s.The  presented  results  clearly  show  that 
the  magnitude  of  the  cell  voltage  overshoot  is  correlated  with 
the  magnitude  of  the  average  current  density  change  (Case  1: 
A V=  0.069  V  and  Case  2:  A V=  0.022  V),  whereas,  the  relaxation 
time  is  almost  independent  of  it. 

In  the  second  strategy,  it  was  assumed  that  the  FU  parameter 
is  constant  during  the  cell  transient  operation.  Therefore,  after  the 
step-load  change  (t= +0),  a  controller  sets  the  fuel  and  air  flow  rates 
proportional  to  that  new  average  current  density,  keeping  the  FU 
and  AER  constant.  Therefore,  in  a  positive  load  change,  both  fuel  and 
air  flow  rates  reach  a  higher  value,  and  vice  versa.  In  Fig.  8  (Case 
3),  an  example  of  the  cell  dynamic  response  under  this  condition 
is  shown.  For  this  example,  the  initial  average  current  density,  fuel 
utilization,  cell  efficiency  and  operating  voltage  are  3000  Am-2  (at 
t  <  0),  FU  =  85%,  77ceii  =  48.1  %  and  V0p  =  0.630  V,  respectively.  At  t  >  0, 
the  average  current  density  is  changed  to  5000 Am-2.  As  can  be 
seen  in  Fig.  8,  there  is  an  undershooting  of  the  cell  voltage  below  the 
new  steady  state  voltage.  However,  the  magnitude  of  this  under¬ 
shoot  is  lower  than  the  overshoot  magnitude,  which  is  observed 
in  Fig.  2,  whereas,  the  variation  of  the  load  change  is  the  same, 
but  in  different  directions.  The  model  output  results  for  V(t  =  +0) 
and  V{t ->  oo),  27Ceii  are  0.50, 0.52  and  39.8%,  respectively.  The  relax¬ 
ation  times  for  recovering  90%,  95%  and  99%  of  the  dynamic  voltage 
drop  respectively  are,  r90%  =  130  s,  r95%  =  155  s  and  r99%  =  190  s.  The 
new  relaxation  data  are  almost  halved,  compared  to  the  previous 
results. 

In  Fig.  9,  the  maximum  and  minimum  temperatures  of  the  cell 
as  well  as  the  maximum  cell  temperature  gradient  for  the  above 
example  are  shown.  In  comparison  to  the  curves  presented  in 
Fig.  6,  both  maximum  allowable  local  temperature  gradient  and 
maximum  cell  temperature  experience  a  lower  change  during  the 
transient  operation  time. 
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Fig.  6.  Maximum  and  minimum  temperatures  and  maximum  cell  temperature  gra-  Fig.  8.  Transient  responses  of  cell  voltage  and  power  to  step-load  change  from 
dient  as  a  function  of  time  (Case  1).  3000Airr2  to  5000  Am-2  (Case  3). 
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Fig.  9.  Maximum  and  minimum  temperatures  and  maximum  cell  temperature  gra¬ 
dient  as  a  function  of  time  (Case  3). 


It  can  be  concluded  that  in  the  second  strategy,  the  cell  response 
is  faster  and  the  presented  strategy  can  be  effectively  used  to  con¬ 
trol  an  SOFC  operation  during  a  load  change. 

Finally,  it  is  worth  noting  that  the  developed  dynamic  model 
shows  the  same  general  response  and  trends  as  the  previous  litera¬ 
ture  models  [6,9,21  ].  The  obtained  relaxation  times  are  also  almost 
similar  to  those  published  for  the  planar  SOFC  cell  [6,9,10,15]. 

4.2.  Transient  response  to  an  inlet  flow  temperature  change 

The  cell  performance  is  highly  dependent  on  its  inlet  flow  tem¬ 
peratures.  As  a  result,  in  a  steady  state  operation  point,  the  fuel 
and  air  temperatures  to  a  stack  are  usually  constant.  However,  in 
some  cases,  i.e.,  due  to  the  failure  of  the  system  components,  such 
as  reformer,  considerable  changes  in  the  inlet  flow  temperatures 
can  occur. 

In  this  section,  the  transient  behaviour  of  a  planar  SOFC  under  an 
inlet  flow  temperature  change  is  investigated  (Case  4).  Therefore,  it 
is  assumed  that  the  inlet  air  and  fuel  temperatures  decrease  about 
50 1<  from  their  initial  values,  i.e.,  1173  K.  In  the  initial  steady  state 
operation,  the  cell  average  current  density,  FU,  AER,  are  3000  A  m-2, 
85%,  7,  respectively.  These  values  are  kept  constant  during  the 
simulation.  The  cell  initial  efficiency  and  voltage,  respectively,  are 
48.1%  and  0.63. 

Fig.  10  shows  the  current  density  distribution  along  the  cell 
length  as  a  function  of  time.  In  comparison  to  Fig.  3,  the  current 
density  variation  is  smoother  here.  The  cell  voltage  and  power 
versus  time  are  shown  in  Fig.  11.  As  can  be  seen,  these  parameters 
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Fig.  10.  Current  density  distribution  along  cell  length  as  a  function  of  time  (Case  4). 
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Fig.  11.  Transient  responses  of  cell  voltage  and  power  to  inlet  flow  temperatures 
change  from  1173  K  to  1123  K  (Case  4). 
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Fig.  12.  Maximum  and  minimum  temperatures  and  maximum  cell  temperature 
gradient  as  a  function  of  time  (Case  4). 


gradually  change  over  time,  without  any  jump.  The  model  out¬ 
put  results  for  V{t^  oo)  and  ?7Cell  are  0.57  and  43.6%,  respectively. 
The  relaxation  time  for  recovering  90%,  95%  and  99%  of  the  volt¬ 
age  drop,  respectively,  are  r90%  =  590  s,  r95%  =  685  s  and  r99%  =  800  s. 
These  new  relaxation  data  are  at  least  two  times  higher,  compared 
to  the  previous  results.  Therefore,  it  can  be  concluded  that,  for  the 
case  investigated,  the  cell  response  time  to  a  flow  inlet  temper¬ 
ature  change  is  slower  than  a  step-load  change.  However,  this  is 
not  generally  correct  as  it  is  dependent  on  several  parameters  i.e., 
magnitude  of  change. 

Maximum  and  minimum  temperatures  of  the  cell,  as  well  as  the 
maximum  cell  temperature  gradient,  are  shown  in  Fig.  12.  In  com¬ 
parison  to  the  previous  results,  there  is  not  any  significant  change 
in  the  maximum  cell  temperatures,  whereas,  there  is  an  increase 
in  the  maximum  cell  temperature  gradient  curve  and  the  cell  tem¬ 
perature  drop  is  more  noticeable. 

5.  Conclusions 

As  an  incremental  step  toward  a  more  comprehensive  under¬ 
standing  of  the  dynamic  behaviour  of  planar  DIR-SOFC  systems,  a 
previously  developed  steady  state  SOFC  cell  model  was  converted 
to  a  dynamic  cell  model.  This  dynamic  model  was  employed  here 
to  explore  the  dynamic  characteristics  of  a  SOFC  cell  under  sev¬ 
eral  step-load  and  an  inlet  flow  temperature  changes.  The  step-load 
changes  were  investigated  under  two  control  strategies.  In  the  first 
one,  it  was  assumed  that  the  fuel  and  air  flow  rates  to  the  cell  were 
constant  during  the  cell  transient  operation.  In  the  second  strat¬ 
egy,  however,  the  fuel  and  air  flow  rates  were  variable  and  their 
values  were  proportional  to  the  cell  new  average  current  density, 
keeping  the  FU  and  AER  constant.  The  results  show  that  the  relax¬ 
ation  time  is  strongly  related  to  the  thermal  behaviour  of  the  cell 
and  the  applied  control  strategy.  However,  it  is  almost  indepen¬ 
dent  of  the  load  variation  magnitude.  It  was  observed  that  the  cell 
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temperature  and  its  gradient  variations,  as  well  as  the  relaxation 
time,  were  smaller  in  cases  where  the  cell  worked  under  a  constant 
FU  during  the  transition  time.  Therefore,  the  second  strategy  can  be 
effectively  used  to  control  an  SOFC  operation  during  a  load  change. 
Results  obtained  from  the  simulation  of  the  cell  transient  behaviour 
during  an  inlet  flow  temperature  change  also  proved  that  for  the 
case  simulated,  relaxation  times  were  at  least  two  times  higher, 
compared  to  the  previous  results  for  the  load  changes. 
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